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Abstract 

(^ ■ We study the valence approximation in lattice QCD of hadrons where the 

cloud quarks and antiquarks are deleted by truncating the backward time prop- 
agation (Z graphs) in the connected insertions. Whereas, the sea quarks are 
eliminated via the quenched approximation and in the disconnected insertions. 
It is shown that the ratios of isovector to isoscalar matrix elements in the nu- 
cleon reproduce the SU(6) quark model predictions in a lattice QCD calculation. 
We also discuss how the hadron masses are affected. 



<D 

1 Introduction 

X 

In addition to its classification scheme, the quark model is, by and large, quite suc- 
cessful in delineating the spectrum and structure of mesons and baryons. One often 
wonders what the nature of the approximation is, especially in view of the advent 
of quantum chromodynamics (QCD). In order to answer this question, we need to 
understand first where the quark model is successful and where it fails. 

To begin with, we need to define what we mean by the quark model. We consider 
the simplest approach which includes the following ingredients: 

• The Fock space is restricted to the valence quarks only. 

• These valence quarks, be them the dressed constituent quarks or the bare 
quarks, are confined in a potential or a bag. To this zeroth order, the hadron 
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wavefunctions involving u,d, and s quarks are classified by the totally symmetric 
wavefunctions in the flavor-spin and orbital space according to the SU(Q) x 0(3) 
and totally antisymmetric/symmetric in the color space for the baryons/mesons. 

• The degeneracy within the the multiplts are lifted by the different quark masses 
and the residual interaction between the quarks which is weak compared to the 
confining potential. The one-gluon exchange potential is usually taken as this 
residual interaction to describe the hyper-fine and fine splittings of the hadron 
masses. 

Given what we mean by the quark model, it is easier to understand where the 
quark model succeeds and fails. It is successful in describing hadron masses, relations 
of coupling and decay constants, magnetic moments, Okubo-Zweig rule, etc. It is 
worthwhile noting that all these are based on the valence picture aided with SU(6) x 
0(3) group for its color-spin and space group. On the other hand, it fails to account 
for the U(l) anomaly (the rj' mass) , the proton spin crisis and the irNa term. All 
these problems involve large contribution from disconnected insertions involving sea- 
quarks [[I]]. It is natural not to expect the valence quark model to work. There are 
other places where the valence quark model does not work well. These include irir, irN 
scatterings, current algebra relations, and the form factors of the nucleon which are 
better described by meson effective theories with chiral symmetry taken into account. 
For example, the im scattering is well described in the chiral perturbation theory, the 
tcN scattering and the nucleon electromagnetic , axial, and pseudoscalar form factors 
(especially the neutron charge radius), Goldberg- Treiman relation are all quite well 
given in the skyrmion approach |2|. One common theme of these models is the chiral 
symmetry which involves meson cloud and hence the higher Fock space beyond the 
valence. 



2 Valence Approximation 

It is then clear that there are three ingredients in the classification the quarks, i.e. 
the valence, the cloud, and the sea quarks. The question is how one defines them 
unambiguously and in a model independent way in QCD. It has been shown recently 
H that in evaluating the hadronic tensor in the deep inelastic scattering, the three 
topological distinct contractions of quark fields lead to the three quark-line skeleton 
diagrams. The self-contraction of the current leading to a quark loop is separated from 
the quark lines joining the nucleon interpolating fields. This disconnected insertion 
(D.I.) refers to the quark lines which are of courses connected by the gloun lines. 
This D.I. defines the sea-parton. One class of the connected insertion (C.I.) involves 
an anti-quark propagating backwards in time between the currents and is defined as 
the cloud anti-quark. Another class of the C.I. involves a quark propagating forward 



in time between the currents and is denned to be the sum of the valence and cloud 
quarks. Thus, in the parton model, the antiquark distribution should be written as 

q\x)=7f c (x) + q i s (x)- (1) 

to denote their respective origins for each flavor i. Similarly, the quark distribution 
is written as 

q\x) = q i v {x) + q i c {x)+q i s {x) (2) 

Since q\{x) = q^x), we define q l c {x) = q l c (x) so that q\r(x) will be responsible for the 
baryon number, i.e. / uy(x)dx = f[u(x) — u(x)]dx = 2 and J dy(x)dx = f[d(x) — 
d(x)] = 1 for the proton. 

We can reveal the role of these quarks in the nucleon matrix elements which 
involve the three-point function with one current. The D.I. in the three-point function 
involves the sea-quark contribution to the m.e. It has been shown that the this 
diagram has indeed large contributions for the flavor-singlet scalar and axial charges 
Q so that the discrepancy between the valence quark model and the experiment in 
the 7rNa term and the flavor-singlet Qa can be understood. Thus we conclude that 
in order to simulate the valence quark model, the first step is to eliminate the quark 
loops. This can be done in the well-known quenched approximation by setting the 
fermion determinant to a constant. 

In order to reveal the effect of the cloud degree of freedom, we have calculated 
the ratios of the isoscalar to isovector axial and scalar charges in a quenched lattice 
calculation. The ratio of the isoscalar (the C.I. part) to isovector axial charge can be 
written as 
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where Au(Ad) is the polarized parton distribution of the u(d) quark and antiquark 
in the C.I. For the non-relativistic case, g\ is 5/3 and g\ for the C.I. is 1 Thus, the 
ratio Ra should be 3/5. Our lattice results based on quenched 16 3 x 24 lattices with 
/3 = 6 for the Wilson k ranging between 0.154 to 0.105 which correspond to strange 
and twice the charm masses are plotted in Fig. 1 as a function of the quark mass 
ma = ln(AK c / ' K — 3). We indeed find this ratio for the heavy quarks (i.e. k > 0.133 
or ma > 0.4 in Fig.l). This is to be expected because the cloud antiquarks which in- 
volves Z-graphs are suppressed for non-relativistic quarks by 0{p/m q ). Interestingly, 
the ratio dips under 3/5 for light quarks. We interpret this as due to the cloud quark 
and antiquark, since in the relativistic valence quark models (i.e. no cloud nor sea 
quarks) the ratio remains to be 3/5. To verify that this is indeed caused by the cloud 
antiquarks from the backward time propagation, we perform the following approxi- 
mation. In the Wilson lattice action, the backward time hopping is prescribed by the 



term — «(1 — ^fA)U^(x)8 Xt y- ai . We shall amputate this term from the quark matrix 
in our calculation of the quark propagators. As a result, the quarks are limited to 
propagating forward in time and there will be no Z-graph and hence no cloud quarks 
and antiquarks. The Fock space is limited to 3 valence quarks. Thus we shall refer to 
this as the valence approximation and we believe it simulates what the naive quark 
model is supposed to describe by design. After making this valence approximation 
for the light quarks with k = 0.148,0.152, and 0.154 (The quark mass turns out to 
differ from before only at the perturbative one-loop order, i.e. 0(a s ), which is very 
small, we find that the ratio Ra becomes 3/5 with errors less than the size of the 
circles in Fig. 1. Since the valence quark model prediction of Ra is well reproduced 
by the valence approximation, we believe this proves our point that the deviation of 
Ra from 3/5 in Fig. 1 is caused by the backward time propagation, i.e. the cloud 
quarks and antiquarks. 

Similar situation happens in the scalar matrix elements. In the parton model 
description of the forward m.e., the ratio of the isovector to isoscalar scalar charge of 
the proton for the C.I. is then approximated according to eqs. (ID) and (0) as 
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Since the quark/ ant iquark number is positive definite, we expect this ratio to be < 
1/3. For heavy quarks where the cloud antiquarks are suppressed, the ratio is indeed 
1/3 (see Fig. 2). For quarks lighter than k = 0.140, we find that the ratio is in fact 
less than 1/3. The lattice results of the valence approximation for the light quarks, 
shown as the circles in Fig. 2, turn out to be 1/3. This shows that the deviation of 
Rs from 1/3 is caused by the cloud quarks and antiquarks. With these findings, we 
obtain an upper-bound for the violation of GSR ||, i.e. n n — ng< —0.12 ±0.05. This 
clearly shows that n s — ng is negative and is quite consistent with the experimental 
result Jdx[u p (x) - d p (x)] = -0.14 ± 0.024. 



3 Hadron Spectroscopy 

To further explore the consequences of the valence approximation, we calculate the 
baryon masses. Plotted in fig. 3 are masses of A, N, p, and 7r as a function of the 
quark mass ma on our lattice with quenched approximation. We see that the hyper- 
fine splittings between the A and N, and the p and ti grow when the quark mass 
approaches the chiral limit as expected. However, it is surprising to learn that in the 
valence approximation, the A and N become degenerate within errors, so do the p and 
tc as shown in Fig. 4. Since the one-gluon exchange is not switched off in the valence 
approximation, the hyper-fine splitting is probably not due to the one-gluon exchange 
potential as commonly believed. Since this is a direct consequence of eliminating the 



cloud quark/ ant iquark degree of freedom, one can speculate that it has something to 
do with the cloud. It seems that a chiral soliton like the skyrmion might delineate a 
more accurate dynamical picture than the one-gluon exchange spin-spin interaction. 

To conclude, we find that the valence approximation in QCD reproduces the SU(6) 
results of the valence quark model better than we anticipated. Especially in hadron 
masses, the results seem to indicate that there are no hyper-fine splittings, modulo 
the uncertainty due to the statistical and systematic errors. 
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Figure Captions: 

Fig. 1 The ratio Ra of eq. (3) as a function of the quark mass ma = Iti{Ak c / K — 3). 

Fig. 2 The ratio Rs of eq. (4) as a function of the quark mass ma. 

Fig. 3 Masses of A, N, p, and 7r (in lattice units) as a function of the quark mass ma 

in the quenched approximation. 

Fig. 4 The same as in Fig. 3 with the valence approximation. 



This figure "figl-l.png" is available in "png" format from: 
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